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Abstract 

A search is performed for heavy long-lived charged particles using 3.0 fb“^ of proton- 
proton collisions collected at ^/s = 7 and 8TeV with the LHCb detector. The 
search is mainly based on the response of the ring imaging Cherenkov detectors to 
distinguish the heavy, slow-moving particles from muons. No evidence is found for 
the production of such long-lived states. The results are expressed as limits on the 
Drell-Yan production of pairs of long-lived particles, with both particles in the LHCb 
pseudorapidity acceptance, 1.8 < r] < 4.9. The mass-dependent cross-section upper 
limits are in the range 2-4 fb (at 95% CL) for masses between 124 and 309GeV/U. 
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1 Introduction 


Several extensions of the Standard Model (SM) propose the existence of charged massive 
stable particles (CMSP). Stable particles, in this context, are long-lived particles that 
can travel through a detector without decaying. These particles can have long lifetimes 
for a variety of reasons, e.g. a new (approximately) conserved quantum number, a weak 
coupling or a limited phase space in any allowed decay. In supersymmetric (SUSY) models, 
CMSPs can be sleptons {£), charginos, or R-hadrons. R-hadrons are colourless states 
combining squarks (^ or gluinos (g) and SM quarks or gluons. In the gauge-mediated 
supersymmetry breaking (GMSB) model Bi the breakdown of SUSY is mediated by 
gauge interactions and can occur at a relatively low energy scale. For a particular range 
of parameter space in the minimal model (mGMSB) the next-to-lightest supersymmetric 
particle can be a long-lived stau (fq), with a mass of the order of lOOGeV/c^ or higher. 
The Ti is the lightest mass eigenstate, resulting from the mixture of right-handed and 
left-handed superpartners of the r, dominated by the right-handed component. 

A GMSP loses energy mainly via ionisation; strongly interacting GMSPs are not 
considered here. In a detector such as LHGb a GMSP with a kinetic energy above about 
5 GeV should be able to traverse the muon chambers. Those particles would often be 
produced with a relatively low velocity and could be identihed by their time-of-flight, and 
by their specihc energy loss, dE/dx, in the detectors; Cherenkov radiation would be absent 
in Cherenkov counters tuned for ultra-relativistic particles. 

Several experiments have searched for CMSPs |4- 12 . With the exception of DELPHI j^, 
which had Cherenkov counters, the analyses are based on dE/dx and time-of-flight 
measurements. The primary interest here is to show the potential of the identihcation 
technique based on ring imaging Cherenkov (RICH) detectors, in addition to the exploration 
of the forward pseudorapidity region only partially covered by the central detectors at the 
Tevatron and the LHC. 

The analysis described in this study is mainly based on the absence of Cherenkov 
radiation in the RICH detectors. This technique is used to search for pairs of CMSPs in 
LHCb, produced by a Drell-Yan mechanism. 


2 The LHCb detector and the detection of slow par¬ 
ticles 


The LHCb detector 13,14 is a single-arm forward spectrometer covering the approximate 


pseudorapidity range 1.8 < rj < 4.9, designed for the study of particles containing 6 or c 
quarks. The detector includes a high-precision tracking system consisting of a silicon-strip 
vertex detector (the vertex locator, VELO) surrounding the proton-proton interaction 
region |^, a large-area silicon-strip detector located upstream of a dipole magnet with a 
bending power of about 4Tm, and three stations of silicon-strip detectors and straw drift 
tubes 16 placed downstream of the magnet. The tracking system provides a measurement 
of momentum, p, of charged particles with a relative uncertainty that varies from 0.5% 
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Figure 1: CMSP velocity spectrum for the CMSP masses of 124GeV/c^ and 309GeV/c^. The 
proton-proton centre-of-mass energy is 7 TeV. The dots with error bars show the efficiency to 
detect tracks as a function of the f3 of the particle (right scale). 


at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary 
vertex (PV), the impact parameter (IP), is measured with a resolution of (15 + 29/pt) 
where pt is the component of the momentum transverse to the beam, in GeV/c. Photons, 
electrons and hadrons are identihed by a calorimeter system consisting of scintillating- 
pad and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. 
Muons are identihed by a system composed of alternating layers of iron and multiwire 
proportional chambers 
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Different types of charged particles are distinguished using information from two RICH 
detectors 18 . The RICH system, which plays a crucial role in this analysis, consists of an 


upstream detector with silica aerogel and C4F10 gas radiators, positioned directly after 
the VELO, and a downstream detector with a CF4 gas radiator, located just after the 
tracking system. 

The online event selection is performed by a trigger 


19 , which consists of a hardware 


stage, based on information from the calorimeter and muon systems, followed by a software 
stage, which applies a full event reconstruction. 

The analysis presented here is based on two data sets collected in 2011 and 2012 
corresponding to integrated luminosities of 1.0 fb~^ and 2.0 fb“^ from proton-proton 
collisions recorded at centre-of-mass energies of 7 and 8 TeV, respectively. 

In the production process considered, CMSPs can have velocities /3 = n/c as low as 0.7, 
and their arrival time at the subdetectors can differ by several nanoseconds with respect 
to lighter particles with /? ~ 1. For illustration, the /3 spectrum is shown in Fig. for 
two values of the CMSP mass, at centre-of-mass energy of 7 TeV. The effects of such 
delayed detection on the efficiencies of the subdetectors are determined from simulation in 
which the timing information is modelled according to dedicated electronic measurements 
and tests in beam. The muon chambers have the largest inefficiency for slow-particle 
reconstruction. The maximal delay for a particle to be accepted by the front-end electronics 
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Table 1: Values of the mGMSB A parameters in the SPS7 scenario used in this study, the 
corresponding masses of the ti, and the cross-section of the pair production at next-to- 
leading order. The last two columns give the detector acceptance A. 


A 

nir 

a 

(fb) 

H(%) 

(TeV) 

(GeV/c 2 ) 

7 TeV 

8 TeV 

7 TeV 

8 TeV 

40 

124 

16.90 ± 0.79 

21.20 ± 0.91 

8.3 

9.5 

50 

154 

7.19 ± 0.38 

9.20 ± 0.46 

6.5 

7.7 

60 

185 

3.44 ± 0.20 

4.50 ± 0.24 

5.2 

6.1 

70 

216 

1.79 ± 0.11 

2.39 ± 0.14 

4.3 

5.0 

80 

247 

1.00 ± 0.07 

1.35 ± 0.08 

3.4 

4.1 

90 

278 

0.57 ± 0.04 

0.80 ± 0.05 

2.8 

3.4 

100 

309 

0.34 ± 0.02 

0.49 ± 0.03 

2.3 

2.9 


is 12.5 ns 17 . In the most downstream mnon chamber, this delay corresponds to the 


arrival of a particle with (5 = 0.83. To be identified as a muon, the charged particle must 
be associated with hits in the last four muon chambers, a requirement that particles with 
/3 < 0.8 fail to meet. The large time-of-flight can also bias the reconstructed position of 
the particle passing through the tracker straw tubes, which accept a maximal drift-time of 
about 35 ns for tracks passing close to the straw radius of 2.5 mm [^. These combined 
effects result in a vanishingly small reconstruction efficiency for particles with (3 < 0.8 but 
an efficiency above 95% li j3 > 0.95, as shown in Fig. [I] 


3 Simulation 


3.1 CMSP signal 

The adopted framework is stau pair production, Ti^Ti~^ in mGMSB via a Drell-Yan 
process. Pairs of CMSPs originating from cascade decays of heavier particles are explicitly 
not considered. In the following the symbol t\ is used when the context is explicitly the 
mGMSB model, while GMSP is kept for the more general context. 

The mGMSB model has six parameters [^|^: the SUSY breaking scale (A), the mass 
scale of the SUSY loop messengers {Mm), the number of messenger supermultiplets (Y 5 ), 
the ratio of the vacuum expectation values of the two neutral Higgs fields (tan/?), the sign 
of the Higgs mass parameter (/r), and the parameter Ggrav, which affects the gravitino 
mass. The Spheno3.0 SUSY spectrum generator 20 is used to compute the masses of 


the t/ as a function of the above six parameters. The SPS7 benchmark scenario is 
used to determine the parameter space, where = 3, tan/? = 15, /i > 0, = 2A, and 

the parameter Gg^av = 4000 are fixed. Variation of A then uniquely determines the f/ mass 
and lifetime, which is of the order of 100 ns. In this study the ti is considered stable. 

The predictions for Tq pair production are based on next-to-leading order (NLO) 
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cross-section calculations by the Prospino 2.1 program 
distribution function (PDF) set 


22 using the CTEQ6.6M parton 
23 . These predictions at y/s = 7 and 8 TeV are presented 
in Table The relative theoretical uncertainties vary between 5% and 8%, and are 


determined following Ref. 24 


Fully simulated signal samples, with masses varying from 124 to 309 GeV/c^, have 
been produced for proton-proton collisions at y/s = 7 and 8 TeV. The t) pairs generated 
by Pythia 6.423 pS], with both Ty particles in the hducial range 1.8 < t] < 4.9 are passed 


to Geant4 1^,27 for detector simulation. The fraction of t) pairs within the hducial 
range is dehned as the acceptance, A. The acceptance factor obtained from Pythia with 
the MSTW2008 PDF set is also shown in Table with model uncertainties ranging 
from 5% to 9% for Ty mass from 124 GeV/c^ to 309 GeV/c^, mainly associated to the choice 
of PDF. 

For larger 7y masses, the Drell-Yan process results in a lower forward boost of the ly 
pair, with a subsequent increase in the pair opening angle in the detector frame. The 
decrease of A for an increasing 7y mass is due to a higher probability for one of the particles 
to escape the LHGb geometrical acceptance. 


3.2 Background 

The main background is from the Drell-Yan production of muon pairs, Z/y* —>■ 
Samples of Z/y* —)■ ii~ events have been produced with Pythia and fully simulated with 
Geant4. The cross-section for this process has been calculated with DYNNLO at 
next-to-next-to-leading order with the MSTW2008 PDF set. The preselection requirements 


(see Sect. 4.1) lead to values of the predicted cross-section in LHGb for y/s =7 and 8 TeV 
of 1.08 ± 0.10 and 1.36 ± 0.12 pb, respectively. These values are nearly two orders of 
magnitude larger than the predicted ly pair cross-section in the most favourable case, 
corresponding to A = 40 TeV. 

Other background sources include muons produced by top quark pairs, and from r pairs. 
To study the background contributions from these processes, samples of Z/y* —)■ 
and top quark pair decays have been simulated. 


r 


4 Data selection 

The event selection is performed in two steps: a preselection aimed at suppressing the 
most prominent backgrounds, followed by a multivariate analysis, based on an artihcial 
neural network that is trained using calibrated simulation. 

4.1 Preselection 

GMSP candidates are identihed as high-momentum charged particles with hits in the 
VELO, all the tracking stations and the four last muon detectors. 

Events are selected that contain two or more such particles where one of the particles 
passes the high-px single muon trigger with a threshold of 15 GeV/c. The trigger efficiency 
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is estimated from simulation to be 92% for a mass of 124GeV/c^, and 89% for 309GeV/c^. 
The two candidates must have opposite charge and each have pt > 50 GeV/c. To reject 
background from Z/ 7 * —)• the pair must have a dimuon mass larger than 100 GeV/c^. 

A mass-dependent lower threshold on momentum is applied to select particles with {3 > 0.8. 

Several criteria are used to reject muons, electrons and hadrons. Pions and kaons 
in jets may be identihed as muons if they decay in flight or if shower fragments escape 
from the calorimeters to the muon stations. As hadrons and electrons deposit more 
energy in the calorimeters than that expected for GMSPs, an efficient rejection of these 
backgrounds is achieved by requiring the sum of the ECAL and HCAL energies associated 
with the extrapolation of the charged particle to the calorimeters to be less than 1 % of 
the momentum of that particle. The background from misidentihed muons contributes 
approximately equally to same- and opposite-charge pairs (^. No same-charge event is 
found in the preselected data, showing that this contribution is negligible. 

GMSPs, as well as muons from ^’/y* decays, would be produced at the PV and should 
have a smaller IP with respect to the PV than muons from heavy quark or tau decays. 
Requiring an IP of less than 50 /xm selects efficiently CMSP candidates. After preselection, 
the contribution from the Z/ 7 * —)■ process where both taus decay leptonically to 

muons is estimated from simulation to contribute less than 0.1 events in total. Pairs 
of muons produced from top quark decays into b quarks and bosons, with the 
bosons decaying leptonically into muons, contribute less than one event, as determined 
from simulation. 

In summary, after preselection the only signihcant source of background is from 
Z/ 7 * —)■ /i^/x“ decays. The predicted number of dimuon events in the 7TeV ( 8 TeV) 
data set is 249±49 (570±110) which is in good agreement with the 239 (713) observed 
candidate events. The uncertainties comprise contributions from the preselection cuts 
(Section]^, and the uncertainty on the ^’/y* —)■ cross-section (Section |3.2[ ). The 

expected number of events with 77 pairs is 2.7 events in the full data set of £ = 3.0fb~^, 
according to the cross-section calculated with Prospino2.1, with SPS7 parameters and a 
77 mass of 124GeV/c^. 

4.2 Selection 

An artificial neural network (ANN) is used to distinguish GMSPs from muons by exploiting 
the difference in interactions that these particles have in matter. To reduce model 
dependence, the ANN is applied to the individual CMSP candidates, rather than to CMSP- 
pairs, and a minimum requirement is placed on the product of the two ANN responses. 
Four variables of the CMSP candidates are used as ANN inputs, computed from the energy 
deposited in the VELO sensors (AE VELO), in the ECAL (AE ECAL), in the HCAL (AE 
HCAL), and a likelihood variable associated with the RICH information (DLLx). Model 
dependence is reduced as much as possible by the absence of kinematical observables 
in the ANN. The energy loss of a charged particle traversing a VELO sensor follows a 
Landau distribution. The most probable energy deposition in a sensor is estimated using 
a truncated mean where only the 60% lowest depositions are averaged. 
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Table 2: Number of events with both CMSPs candidates satisfying DLLx > — 5, and the hnal 
efficiency, e, after the multivariate analysis selection, given for each mass hypothesis. 


^CMSP 
(GeV/c^) 

CMSP-pair 

candidates 

7 TeV 8 TeV 

e 

7 TeV 

(%) 

8 TeV 

124 

38 

73 

49.6 ± 4.4 

45.1 ± 4.4 

154 

36 

68 

48.9 ± 4.5 

44.5 ± 4.5 

185 

36 

68 

46.0 ± 4.7 

41.9 ± 4.7 

216 

28 

56 

42.0 ± 4.8 

38.5 ± 4.8 

247 

24 

49 

37.5 ± 5.0 

35.0 ± 5.0 

278 

24 

49 

32.8 ± 5.1 

31.2 ± 5.1 

309 

13 

30 

28.4 ± 5.3 

27.3 ± 5.3 


Particle identification for a CMSP candidate, using RICH information, is provided 
by the DLLx variable. A particle identification hypothesis is assigned to a track using a 
likelihood method. The information from the three radiators is combined and a “delta 
log-likelihood” (DLL) value computed. The DLL gives, for each track, the change in the 
overall event log-likelihood when the particle ID hypothesis is changed from vr to fi, e, 
K, p. The DLLx classihcation has been added to account for high momentum particles 
which do not radiate, or have a Cherenkov angle which is too small to £t one of the hve 
particle hypotheses. A positive DLLx indicates a high probability that the candidate has 
a relatively low velocity. More details are given in Section 

Events with both candidate CMSPs with DLLx > —5 are used in the analysis, with 
no loss of signal, as deduced from simulation. The numbers of selected events with 
CMSP-pairs are given in Table 

Simulated events are used to train the ANN. The hrst three variables dehned above are 
calibrated using muons from Z/y* —>■ to ensure that simulation agrees with data. A 

total of 25k Z events from the 2011 data set and 65k from the 2012 data set are used for 
the calibration. In the Z mass region the expected amount of signal is smaller than one 
event and cannot bias the procedure. The DLLx variable is by far the most discriminating, 
and its calibration procedure is presented in detail in Section 

The ANN training is carried out independently for the 7 and 8 TeV data sets, and 
for all the CMSP mass hypotheses considered. Figure shows the distribution of the 
four ANN input variables for the 8 TeV data set, compared to the background and signal 
predictions; good agreement can be seen between data and simulated background. 

The discriminating variable is the product of the ANN outputs obtained for the two 
CMSP candidates. This “pair signihcance” is shown in Fig. for the 124 and 309GeV/c^ 
CMSP mass hypotheses for both 7 and 8 TeV data sets. The requirement placed on the 
pair signihcance is determined by the value needed to achieve a 95% signal efficiency. 
After the pair signihcance selection, the signal efficiency for candidate events in the LHCb 
acceptance is 50% for CMSPs with a mass of 124 GeV/c^, decreasing for increasing mass to 
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Figure 2: Number of CMSP candidates, as a function of the four variables used as inputs to 
the ANN. There are two CMSP candidates per event. The black dots with error bars show the 
2012 data. The dashed red histogram is the expected shape for 124GeV/c^ CMSPs and the blue 
histogram shows the background from Zj^* decays into muons. The energy in the VETO is 
given in units of minimum ionising particle (MIP) deposition. The first bin of the histogram for 
AE in the ECAL has been multiplied by a factor 0.25. 


a minimum of 27% at 309 GeV/c^. The signal efficiency values for CMSPs in the acceptance, 
after the ANN selection, are given in Table After the full selection is applied, the 
diniuon background is suppressed by a factor of 10“^. 

5 CMSP identification with Cherenkov detectors 

The present study uses the Cherenkov radiation produced in the RICH detectors to identify 
CMSPs. The Cherenkov momentum thresholds for muons, protons, and CMSPs with 
masses of 124GeV/c^ and 309GeV/c^, are given in Table for the three radiators in the 
LHCb detectors. Only CMSP candidates with momenta above 200 GeV/c are considered. 
For this momentum range, particles with masses of the order of MeV/c^ to GeV/c^, 
have Cherenkov angles very close to the saturation value arccos(l/(n/3)), where n is the 
refractive index of the medium. The fraction of CMSPs with momentum above 2TeV/c is 
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Figure 3: Number of CMSP pairs, N, as a function of the pair significance. The left and right 
figures correspond to the CMSP masses of 124 and 309GeV/c^, respectively. The black points 
with their statistical uncertainty show the 7 TeV (top) and 8 TeV (bottom) data sets. The red 
dashed histogram is the expected shape from CMSP pairs and the blue histogram the background; 
both are normalised to the number of events. The arrows indicate the chosen selection criteria. 


negligible, and the CMSPs are therefore expected not to produce Cherenkov radiation in 
the gaseous radiators. Around half of the 124GeV/c^ CMSPs have a momentum above the 
Cherenkov threshold for aerogel, and Cherenkov angles smaller than the saturation value. 
This allows them to be separated from the background. Only a few percent are expected 
to be in the momentum range 1.4 — 2.0TeV/c, corresponding to Cherenkov angles from 
0.225 to 0.234 rad. It is possible to distinguish these angles from the saturation value of 
0.242 rad in the aerogel as the angular resolution is about 5.6mrad. 

As previously said, the variable DLLx has been introduced to identify high momentum 
particles which do not radiate, or have a Cherenkov angle which is too small to ht one of 
the five particle hypotheses, tt, /i, e, K, p. The DLLx value is positive for the momentum 
distributions of the CMSPs, for all of the masses considered]^ 

^An anomalous signature in the RICH detectors could also be produced by tracks with wrongly assigned 
momenta. This can happen if the particle has an absolute electric charge that is different from the proton 
charge. For instance, a IGeV/c^ proton-like particle is expected to produce Cherenkov light in the two 


8 








































































Table 3: Refractive indices and Cherenkov /3 thresholds for the three radiators. The momentum 
threshold is given for muons, protons, and 124 and 309GeV/c^ CMSPs. 


Radiator 

n 

/^thresh 



Pthresh ( GeV/c) 

h 

P 

CMSP(124) 

CMSP(309) 

Aerogel 

1.03 

0.9709 

0.428 

3.8 

502 

1252 

C 4 F 10 

1.0014 

0.9985 

2.00 

17.7 

2342 

5069 

CF 4 

1.0005 

0.9995 

3.34 

29.7 

3921 

9767 


Simulated dimuon background events and CMSP signal samples used to train the ANN 
are first validated with data. 

The study of the background samples is performed on a set of mnons above the 
Cherenkov threshold and selected from Z decays. Such events have an event topology and 
kinematics that are very close to those of the dimuon background expected in the CMSP 
analysis. The DLLx distribution is shown in Fig. (a) for muons from data and simulated 
Z decays. For illustration, the expected signal shapes for CMSPs with masses of 124 and 
309GeV/c^ are snperimposed. The small difference is dne to a change in the nnderlying 
event and some light from the aerogel for the 124GeV/c^ case. A clear separation between 
the signal and backgronnd mnon DLLx distribntions is observed. The difference in the 
data and simulated muon distribntions is mainly dne to the lack of precision in the 
mapping of the photon detection efficiency in the RIGH system. In particular, the peak 
at DLLx > —5 is produced by the decrease of the photon detection efficiency when 
approaching boundaries in the RIGH modules. The simulation only partially reproduces 
this behaviour and the number of candidates above DLLx = —5 is too low by aronnd a 
factor of two. To compensate for this, 15% of simulated muon events with DLLx falling 
close to zero have been shifted by an ad hoc value to obtain the best agreement between 
data and simnlation in the DLLx > —5 region, resnlting in the distribntion shown in 
Fig. 1^ (b). It is expected that a correct efficiency map should produce such a shift, moving 
above zero the slightly negative DLLx values. This set of simulated background events is 
used to train the ANN. Note that only candidates with DLLx > — 5 are used in the ANN. 
In order to assess the systematic uncertainty associated to this correction method, two 
other procedures are considered. In the hrst procedure, the ANN training is performed on 
the original background simulated data set. In the second, the DLLx values for each muon 
are randomly chosen following a set of templates inferred from the DLLx data distributions 
as fnnction of p and p. Despite the fact that this operation is done in bins of p and p, it 
is obvions that most of the correlation is lost in the randomisation process. The three 
methods are found to provide the same hnal discrimination power and their contributions 
to the systematic nncertainties are small. This is dne to the strong separation between 

gaseous radiators when the measured momentum is above 30GeV/c. If the particle has one third of the 
proton charge, the measured momentum is overestimated by a factor of three and this will lead to an 
incorrect calculation of the Cherenkov emission. 
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Figure 4: (a) DLLx distributions from muons selected from Z decays in 7 TeV data (black 
points) and simulation (blue histogram). The expected shape from 124GeV/c^ and 309GeV/c^ 
CMSPs are also shown, (b) The simulated DLLx values are shifted in such a way that data and 
simulation have the same fraction of entries with DLLx > — 5. 


signal and backgronnd that is provided intrinsically by the DLLx variable. 

The validation of the signal sample is more complex due to the absence of a SM 
process that can be used for calibration. The quality of the simulation was studied 
using protons from T —pvr decays with a velocity below the Cherenkov threshold. The 
differences between the data and simulation for these protons are extrapolated to the 
CMSP kinematics using a fast simulation method, and the contribution to the systematic 
uncertainty estimated. 

The proton is below the Cherenkov threshold in all of the RICH radiators for p < 
3.8GeV/c, and above the threshold for p > 30GeV/c. Pairs of opposite-charge tracks are 
selected from data and from simulated events passing a minimum bias trigger. The pair 
must combine to form a particle with a mass compatible with the known mass of the A 
baryon, and the reconstructed vertex must be more than 3 mm from the beam axis. 

Samples of protons below and above Cherenkov threshold are obtained by choosing the 
momentum regions below 3.8GeV/c and above 30GeV/c, respectively. Figure shows the 
corresponding DLLx distributions, indicating a reasonable agreement between the DLLx 
distributions from data collected at 7 TeV and simulation. A Gaussian function plus a 
polynomial, to account for the tail at low DLLx, is htted to the DLLx distributions for 
below-Gherenkov-threshold particles. The mean and width of the Gaussian functions are 
6.5 ± 0.3 and 3.4 ± 0.3 for data, and 6.0 ± 0.3 and 3.6 ± 0.4 for simulation. The DLLx 
value in data is 0.5 ± 0.4 units higher, which may indicate that there is a lower photon 
detection efficiency in data compared to simulation. A maximal deviation of ±1 DLLx 
units is considered in the following to assess the systematic effects. The extrapolation from 
the low momentum proton result to the GMSP regime is made using a fast simulation. 

In addition to the GEANT4-based full simulation, a fast simulation describing the 
main features of the RIGH measurement process is also used. This allows the impact of 
varying parameters and the algorithms to be studied in a more efficient way. The fast 
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Figure 5: The DLLx variable for: (a) protons selected above the Cherenkov threshold, p > 
30GeV/c, in 7 TeV data (black line) and simulation (blue points); (b) and (c) below-threshold 
protons, p < 3.8GeV/c, from data and simulation, respectively. The results from the fit of a 
Gaussian function plus polynomial to the below-threshold proton distributions are shown by the 
red curves. 


simulation generates a target particle (a proton from A decays, a muon or a CMSP) with 
a momentum distribution representing the phenomenon under study. The underlying 
event is represented by a number of pions with a momentum distribution obtained from 
minimum bias events. The simulation of Cherenkov emission in the radiators is then 
performed for each particle. The number of Cherenkov photons generated by a particle of 
velocity (3 follows a Poisson distribution of average No(/9^u^ — l)/(/3^(?7,^ — 1)), where Nq 
is the maximal number of photons for a saturated ring and n is the refractive index. The 
ring has an average radius corresponding to the expected Cherenkov angle and a Gaussian 
profile of width Uc representing the angular resolution of the detector. Finally, random 
noise is added using the probability for a pixel to fire, probnoise- The nominal values of the 
parameters used in the fast simulation are given in Table The event log-likelihood for 
each target particle hypothesis is 

all pixels active pixels 

LL = - i/i+ ^ In(e^'-l) (1) 

pixel! pixel! 
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Table 4: Nominal values of the parameters used in the fast simulation. 


parameter 

aerogel 

C 4 F 10 

CF 4 

n 

1.03 

1.0014 

1.0005 

No 

8 

28 

24 

(Tc [mrad] 

5.6 

1.6 

0.7 

probnoise 

3% 

3% 

3% 


60 


50 


40 


i 30 

c 

LU 

20 


10 


-20 0 20 
DLLx 

Figure 6: DLLx for protons with p < 3.8GeV/c. In black the result with nominal simulation 
parameters are shown, (a) The red, dashed, and green, dash-dot, plots are for a change by 
+20% and -20% of the photon yield; (b) is for a change by -|-40% (red, dashed) and -40% (green, 
dash-dot) of the random noise probability. 




where z/j is the probability to have photons in the pixel i, including the random background. 
Note that the formula is valid for the binary readout implemented in the RICH electronics. 
The centre of each Cherenkov ring is dehned by the true particle direction. The DLL 
values are subsequently computed. 

Simulated distributions for protons from A decays with p < 3.8GeV/c are shown in 
Fig-i The average DLLx is 6.2, with an RMS of 4.0, for a fast simulation made using the 
nominal parameters. Figure (a) also shows the distributions after varying the detection 
efficiency by ±20%. The corresponding distributions are shifted by ±1.2 units. The 
random noise probability was changed by ±40% from its nominal value which produces 
negligible variation as seen in Fig. (b). This study shows that a variation of ±1 DLLx 
units is obtained by changing the photon detection efficiency by ±15%. A variation of the 
same size can be obtained by changing the angular resolution cTc by 50%. 

The DLLx distributions for CMSPs obtained from the fast simulation with nominal 
parameters are consistent with those obtained from full simulation. Changing the photon 
detection efficiency by ±15%, as inferred from the A study, a variation of ±2 DLLx units 
is obtained. An identical result can be obtained by changing the angular resolution. 
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Table 5: Systematic uncertainties (in %) for the selection of the signal of CMSP pairs and on 
the background retention. Where relevant, the lower value corresponds to the CMSP mass of 
124GeV/c^ and the higher one to 309GeV/c^. 



7TeV 

8 TeV 

Signal 

efficiency 

Background 

retention 

Signal 

efficiency 

Background 

retention 

Luminosity 

1.7 

1.7 

1.2 

1.2 

Trigger, reconstruction 

2.0 

2.0 

2.0 

2.0 

Delayed signals 

2.2-3.7 

0 

2.2-3.7 

0 

IP calibration 

0 

0.7 

0 

0.7 

Hadron electron rejection 

0 

15 

0 

15 

Neural Network 

2.9 

50.0 

2.9 

50.0 

Total syst. uncertainty 

4.5-5.4 

52.3 

4.3-5.2 

52.3 


In summary, an uncertainty of two DLLx units is inferred from the comparison of 
below-threshold protons in data and simulation when using A decays. The extrapolation 
to the CMSP regime is obtained by varying the simulation parameters and leads to an 
uncertainty of four DLLx units on the average DLLx value for such particles. 


6 Uncertainties and results 


After the ANN selection the signal prediction for the chosen model is 2.5, 0.9 and 0.3 
events for the 124, 154 and 185GeV/c^ ti masses, and below 0.1 for the other mass values. 
The expected background is negligible, less than 0.02 events for all the mass hypotheses. 

A summary of the systematic uncertainties is given in Table The total systematic 
uncertainties are approximately 5% for the signal yield and 50% for the background yield. 

Two methods are used to determine the luminosity; a Van der Meer scan and a 
beam-gas imaging method 31 . The uncertainties on the integrated luminosities are 1.7% 


for the 7TeV data set and 1.2% for the 8TeV data set. 

The efficiency for triggering, reconstructing and identifying high-pT muons has been 
studied in detail for the LHCb Z and W boson cross-section measurements 32,33 , and the 
agreement between data and simulation was found to be better than 2%. This percentage 
is taken as the corresponding uncertainty for this analysis. 

A further efficiency uncertainty is considered due to the delayed signals in the tracking 
and muon systems. The timing precision affects the amplitude recorded by the front-end 
electronic boards and the measurement of the drift time in the straw tubes. The effect on 
the efficiency due to a timing uncertainty of ±1 ns is determined from simulation as a 
function of the (3 of the particle. Subsequently, a weighted average of the uncertainty is 
obtained from the (3 distributions for each mass hypothesis, providing values varying from 
2.2% to 3.7% for CMSP masses from 124GeV/c^ to 309GeV/c^. 
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The comparison of the IP distributions in data and simulated events 

indicates a maximal discrepancy of ±5 /im. By changing the requirement on the IP param¬ 
eter by this amount, the corresponding efficiency variation is ±0.7% for the background 
and negligible for the signal. 

The hadron and electron rejection is affected by the calibration of the calorimeters. 
From the comparison of Z/7* —)■ decays in data and simulation a relative uncertainty 

of 10% is inferred. This translates into a 15% change on the background yield, while the 
signal is almost unaffected. 

The training of the ANN is affected by the uncertainty on the background and signal 
models. A 2.7% contribution to the signal efficiency uncertainty is associated with 
the calibration procedures, determined by the comparison of data and simulation for 
Z/ 7 * —)■ and A —)■ pTT decays. Error propagation is performed by modifying the 

ANN training sets, while keeping the test sets and the pair signihcance selection hxed. 
Adding the uncertainties in quadrature with the statistical uncertainty of 1%, gives a total 
of 2.9%. 


The ANN selection leaves a very small amount of simulated background. The binomial 
uncertainty on the background retention is large, at approximately 50%. This value is 
assigned as the uncertainty on the background selection efficiency. 

As already stated, the acceptance A is affected by model uncertainties in the range 
from 5% to 9% for 77 mass from 124GeV/c^ to 309GeV/c^. In addition, the choice of the 
PDF affects the efficiency by modifying the momentum of the CMSP. By scanning various 
PDFs, we have found that this effect is small, not larger that 0.4%, for all the models. 

The cross-section upper limits are computed using the Feldman-Gousins method 
for zero observed candidates, taking into account the expected number of background 
events and the uncertainties 35 . The predicted amount of background is so small that 


it has no sizeable effect on the result. The upper limits at a 95% conhdence level (CL) 
for CMSP pair production in the LHCb geometrical acceptance at -y/s = 7 and 8 TeV are 
listed in Tableand shown in Fig. [^together with the theoretical cross-sections calculated 
for the particular model described in Section |3.1[ 


7 Conclusions 

A search for pairs of long-lived charged particles, with masses in the range 124-309 GeV/c^, 
using 77 pairs predicted by the mGMSB model as a benchmark scenario, is performed 
using data from proton-proton collisions at 7 and 8 TeV, corresponding to an integrated 
luminosity of 3.0 fb“^, collected with the LHCb detector in the forward pseudorapidity 
range 1.8 < rj < 4.9. The candidates are assumed to interact only through the electroweak 
interaction in the detector. Hence, they behave like heavy muons and their main signature 
is the absence of a signal in the RICH detectors. The detection efficiency is limited to 
particles with f3 > 0.8 due to the acceptance in time after beam crossing. The main 
background contribution comes from Z/ 7 * —)■ and is reduced to less than ~ 0.02 

events. No events have been observed. Upper limits are set on the Drell-Yan CMSP pair 
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Table 6: Cross-section upper limits at 95% CL for CMSP pair production in the LHCb acceptance 
in the 7 and 8 TeV. 


^CMSP 

(GeV/c2) 

Upper limit (fb) 

7 TeV 

8 TeV 

124 

6.1 

3.4 

154 

6.2 

3.5 

185 

6.6 

3.7 

216 

7.2 

4.0 

247 

8.1 

4.4 

278 

9.2 

5.0 

309 

10.7 

5.7 



CMSP mass (GeV/c^) 

Figure 7: Upper limits at 95% CL on the cross-sections for the pair production of CMSPs in the 
LHCb acceptance (points) and the corresponding predictions assuming the Drell-Yan production 
of Ti (bands representing Tier uncertainty) with SPS7 parameters, for proton-proton collisions 
as a function of the CMSP mass at -y/s =7 and 8 TeV. 


production cross-section. For proton-proton collisions at y/s = 7 TeV, the 95% CL upper 
limits for the production cross-section of a pair of CMSPs in the LHCb acceptance vary 
from 6.f fb for a mass of 124GeV/c^ up to 10.7 fb for a mass of 309GeV/c^. At ^/s = 
8 TeV, they vary from 3.4 fb to 5.7 fb for the same masses. 

In LHGb the identiheation of CMSPs relies on the energy deposited in the subdetectors, 
the main discrimination power being provided by the RICH system. Together with the 
forward pseudorapidity coverage, this unique feature allows LHCb to complement the 
searches undertaken by the central detectors at the Tevatron and LHC. 
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